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Abstract: Supersymmetric monojets may be produced at the Large Hadron Col- 
lider by the process qg —> qx® — > qXxXHt leading to a jet recoiling against missing 
transverse momentum. We discuss the feasibility and utility of the supersymmetric 
monojet signal. In particular, we examine the possible precision with which one can 
ascertain the Xi<7<7 coupling via the rate for monojet events. Such a coupling contains 
information on the composition of the Xi an d helps bound dark matter direct detec- 
tion cross- sections and the dark matter relic density of the Xv It also provides a check 
of the supersymmetric relation between gauge couplings and gaugino-quark-squark 
couplings. 
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1. Introduction 

With the first Large Hadron Collider (LHC) experimental run now underway, searches 
for physics beyond the Standard Model will begin to test new theories of TeV-scale 
phenomena in parameter regimes that were not accessible to previous collider fa- 
cilities. There are many proposals for theories of TeV-scale physics beyond the 
Standard Model. These include supersymmetry, technicolor, little Higgs models, 
extra-dimensions, low-scale quantum gravity, etc., all of which posit the existence of 
new classes of fundamental particles [1]. In many such models, the lightest member 
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of the new class is absolutely stable due to the presence of a conserved discrete sym- 
metry. As a result, production of the lightest new physics particle (either directly 
or at the end of chain of heavier decaying particles) will lead to missing transverse 
energy signatures at the LHC. More generally, signatures from different models of 
TeV-scale physics often possess similar features, and distinguishing among model 
interpretations may be quite challenging. 

In this paper, we assume that after a number of years of LHC running, a class of 
new particles has been discovered. We assume that their masses and Standard Model 
quantum numbers will be relatively well determined. We also assume that we will 
have some information about the spins of the new particles. Under these assumptions, 
how strong will the case be for a supersymmetric (SUSY) interpretation of the new 
physics? 

Ideally, one would first try to identify the new particles as superpartners of 
Standard Model particles. Such an attempt will likely be incomplete, as the entire 
spectrum of the new physics may not have been revealed, and the spins of the new 
particles may not be reliably known in all cases. In isolation, the discovery of a 
neutral color octet fermion or a color triplet scalar does not necessarily imply the 
discovery of the gluino and the squark. Due to supersymmetry-breaking, the masses 
of the new particles would not provide any evidence for a supersymmetric interpre- 
tation. However, supersymmetry-breaking effects will typically have a small impact 
on dimension-four couplings [2]. These couplings therefore reflect the underlying 
supersymmetric structure. Hence, the measurement of a relation among couplings 
could provide very strong evidence for a supersymmetric interpretation of the new 
physics. 

In the supersymmetric extension of the Standard Model, the couplings of glu- 
ons to squark pairs and gluino pairs and the couplings of the photon, W ± and Z 
to squark pairs and slepton pairs are governed by SU(3) xSU(2) xU(l) gauge invari- 
ance. In contrast, supersymmetry relates the gauge couplings of particles (or their 
supersymmetric partners) to the Yukawa couplings of gauginos to particle-sparticle 
pairs. Similarly, supersymmetry relates the Yukawa couplings of Higgs bosons to 
the Yukawa couplings of higgsinos to particle-sparticle pairs. Thus, the precision 
measurement of the interaction strengths of particle-sparticle vertices can provide 
the smoking gun for a supersymmetric interpretation of the new physics. 

Precision measurements of couplings is usually the domain of e + e~ colliders. 
Indeed, there is a rich program of precision supersymmetry proposed for an exper- 
imental program at the ILC [3]. However, until the ILC becomes a reality (current 
projections suggest that this is unlikely during the present decade), we must rely on 
the LHC running at high luminosity to provide the necessary data for interpreting 
the fundamental nature of new physics discoveries. Ultimately, a combined LHC/ILC 
analysis would yield the most precise tests for the supersymmetric interpretation of 
new physics [4]. 
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Until now, only a few analyses have been proposed for the LHC In Ref. [5], 
same flavor squark pair production at the LHC was investigated in order to measure 
the strength of the gluino- squark- quark coupling. But this analysis still relied on 
additional data from a linear collider. The analysis was refined in Ref. [6] by taking 
into account more production channels of strongly interacting SUSY particles. It was 
found that a test of the SUSY relation for the gluino-squark-quark operator based 
solely on LHC data might be possible with a precision better than 10%. A later 
study [7] assuming 300 fb _1 of LHC data, with varying degrees of ILC input, found 
similar fractional precisions of 3-7%. Finally, in Ref. [8], (left-handed) squark pair 
production was considered via color-singlet gaugino exchange in the t- and u-channel. 
It was shown that a measurement of this process, and therefore a measurement of 
the wino-squark-quark coupling, might be possible at the LHC. 

We propose the development of a program of precision measurements of new 
particle interactions at the LHC in order to provide definitive evidence in support of 
a possible supersymmetric interpretation of TeV-scale physics. Such a program must 
necessarily comprise a broad class of new particle signatures. Ultimately, a global fit 
to a plethora of supersymmetric observables will be required to provide the maximal 
coverage of the underlying parameter space. In this paper, we take the first step 
by identifying a particular signature that is sensitive to the squark-quark-gaugino 
coupling. 

In Section 0, we identify the supersymmetric monojet signal as a promising 
arena for measuring the squark-quark-gaugino coupling. Two basic supersymmetric 
scenarios are considered, which depend on the identity of the lightest supersymmetric 
particle (LSP). We contrast the phenomenology of a bino-LSP and a wino-LSP, which 
arise in many models of TeV-scale supersymmetry. In Section |||, we outline the 
process by which one can test the supersymmetric coupling relations at the LHC. 
These procedures are then employed in Sections |4] and [5] to measure the squark-quark- 
gaugino couplings in the case of the bino-LSP and wino-LSP, respectively. Finally, 
we state our conclusions in Section ||. The pr-spectrum of the monojet signal exhibits 
a peak at a characteristic value of the transverse momentum that is approximately 
given by 



Jacobian peak of the electron transverse momentum in W production and decay. 





This result, which is derived in Appendix A [cf. eq. ( [A . 6 1| ) ] , is analogous to the 
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2. Supersymmetric Monojets Revisited 



2.1 Associated Production of Squarks and Neutralinos 

Supersymmetric monojet signatures have an interesting history. In the early 1980s, 
CERN pp collisions resulting in single jets ("monojet") in association with miss- 
ing transverse energy in the UAl experiment [9] were famously interpreted in terms 
of supersymmetric particle production [10-22]. One example of such a process was 
squark photino production, followed by decay of the squark into a quark and photino. 
Later, it was calculated that the monojet events would predict more di-jet missing 
transverse momentum events than were observed [23,24]. Finally, backgrounds (par- 
ticularly W production, where W — > tv t , followed by r — > u T j), and Zj production, 
where Z — > vv [25], were seen to adequately account for the monojet events. 

Several authors have discussed the ability of LHC experiments to discover large 
extra dimensions through the monojets signal [26,27] predicted by the production of 
a hard gluon and a Kaluza-Klein tower of gravitons. Obvious Standard Model (SM) 
backgrounds in the LHC environment include di-jet production, where one of the 
jets is lost, and the Wj, Zj backgrounds mentioned above. We wish to measure the 
rate of the dominant SM background (i.e. monojet production via Zj) by measuring 
the cases where Z — > (£ = e, /i) and using them to predict the background 

process. Isolated lepton vetoes help discriminate against Wj, and a large missing 
Pt cut discriminates against the QCD di-jet background. More recently, unparticle- 
jet production, leading to a monojet signature has also been examined [28]. The 
monojet pt spectrum is fairly featureless in either the large extra dimensional case 
or in the unparticle-jet case, although the rates are in general different. 

Early -y/i = 7 TeV LHC collisions resulting in monojets recorded by ATLAS [29] 
have already been made public. The data show no excess over SM Monte Carlo 
predictions for 70 nb -1 of LHC data, but we find that the monojet channel is not 
yet sensitive to SUSY models that are not already ruled out by other experiments. 1 
Much more data will be required. 

We wish to appraise the ability of the LHC to measure SUSY monojet produc- 
tion, assuming the minimal supersymmetric standard model (MSSM) via the process 
qg — > qxi- At tree-level, there are two diagrams that contribute, which are exhib- 
ited in Fig. [I]. The tree-level scattering amplitudes for qg — > qxi were obtained in 
Ref. [31]. Strong SUSY-QCD corrections and the leading logarithmic electroweak 
corrections at one loop have been treated in Ref. [32] (although we note that the 
one- loop QCD corrections to qg — > qxi do not yet appear in the literature). 

The processes depicted in Fig. [I] do not constitute the best SUSY search chan- 
nel, since the production amplitude is proportional to a weak (A) times a strong 

1 Note, however, that multijets plus jj/ T channels are already becoming competitive with published 
Tevatron bounds in certain corners of MSSM parameter space [30]. 
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Figure 1: Tree- level Feynman diagrams leading to qg — > qx\- Monojet signatures result 
from q — )■ qx\- 

gauge coupling. Di-squark, di-gluino and/or gluino-squark production are expected 
to have much higher rates because their amplitudes are proportional to the strong 
gauge coupling squared. This is probably why the associated production of squarks 
and neutralinos has not been extensively examined in the literature. On the other 
hand, if the processes in Fig. |l] could be identified and the rates measured, an es- 
timate for the coupling A may result. This coupling contains information on the 
identity of the lightest neutralino x? ; which may be dominantly bino, wino, higgsino 
or a mixture. For example, in the bino (wino) dominated case, supersymmetry dic- 
tates that A is proportional to g' (g), the gauge coupling of the gauge group U(l)y 
(SU(2)i), respectively (where the proportionality constant is fixed by supersymme- 
try). If one knew the constitution of the neutralino from other measurements, then 
the measurement of the coupling A in associated squark-neutralino production would 
provide a direct test of supersymmetry. 2 

It has long been postulated that Xi constitutes the dark matter of our universe. 
Knowing its couplings is a vital ingredient in calculating how much of it is left as a 
thermal relic in the universe. For example, neutralinos may annihilate in the early 
universe through the process depicted in Fig. ||a before freezing out. Such a process 
involves the same vertex as the one in our signal. If Xi constitutes the dark matter 
of the universe, direct detection experiments have the chance to measure it through 
it causing nuclear recoils. As depicted in Fig. ||b, nuclear recoil mediated via squark 
exchange could contribute a significant part of the direct detection cross-section [33] . 
The measurement of A would bound the direct detection cross-section contribution 
from such a channel. 

There are advantages and disadvantages to the SUSY monojet signature as com- 
pared to monojets resulting from large extra dimensions or unparticle-jet production. 
As we shall illustrate below, SUSY monojet production has features in the px spec- 
trum, which would help to convince us that the distribution is not merely an incom- 
plete understanding of SM backgrounds. In particular, the location of the Jacobian 

2 As emphasized in Section 1, the measurement of the gluon couplings to gluino pairs or to squark 
pairs does not constitute a test of supersymmetry, as these latter couplings are governed by QCD. 
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Figure 2: Tree-level Feynman diagrams contributing to (a) dark matter annihilation and 
(b) dark matter direct detection. 

peak provides an indication of the masses of the squark and neutralino (an analytical 
approximation is provided in Appendix A). The disadvantage is that SUSY back- 
grounds may be problematic. Other SUSY processes than those in Fig. p] leading to 
the monojet signature include: X1X1 production, with an initial-state radiated jet, 
and qq production, where each q — > qxl, but both quarks are either in the same 
direction, or one is lost (for example because it is emitted at high rapidity). The 
case where both quarks are in the same direction can be discriminated against with 
a maximum jet mass cut. 

Various papers have examined SUSY monojets recently, but these are all differ- 
ent processes from the one we attempt to isolate and would be classified as SUSY 
backgrounds by our study. In Ref. [34], a scenario where gluinos and neutralinos 
are quasi mass-degenerate were considered. Such a scenario allows a gravitino dark 
matter candidate while being compatible with high reheating temperatures required 
by thermal leptogenesis. This scenario leads to an effective monojet signature via qg 
production, where the gluinos decay to soft QCD radiation and quasi-stable neutrali- 
nos, and the squark decays to a jet and a neutralino. Weak gaugino pair production 
plus a jet, leading to monojet signatures, has recently been examined in Ref. [35] for 
the case of quasi-degenerate gauginos, so that the visible products of their decays are 
too soft to be detected. In Ref. [36], SUSY monojet signatures were examined for a 
region of the MSSM consistent with baryogenesis and dark matter constraints. These 
monojet signatures originated from stop pair production plus an additional QCD jet, 
where the stops decay invisibly into a soft charm quark and a Xi- Monojet searches 
at the Tevatron have recently been used to place bounds on dark matter direct de- 
tection rates [37] and on indirect dark matter searches via gamma ray lines [38]. The 
collider process investigated in Ref. [37] consisted of initial state radiation in dark 
matter pair production, which is classified as a background in the present paper. 
It was found that for very light dark matter, below 5 GeV, the inferred Tevatron 
bounds are stronger than those from direct detection. 3 Pair production of dark mat- 
ter particles in association with one or more jets at Tevatron and LHC, was also 

3 Although the absence of direct detection signals is always subject to a potentially large system- 
atic of unknown density of dark matter at the site of the experiment. 
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considered in Ref. [39]. For our study, all processes mentioned above in this para- 
graph are classified as SUSY backgrounds. An earlier work [40] included predictions 
for monojet signatures from ultra-light gravitinos at the Tevatron and LHC. This is 
a different scenario from the one of interest in the present paper. 

LHC and Tevatron cross sections for total SUSY monojet production were listed 
in Ref. [41], which was primarily concerned with the case of massless Xi particles. 
Of course, the total SUSY monojet production includes our signal process (along- 
side SUSY backgrounds). We extend this study in several ways: for example, we 
perform a more detailed analysis, finding reasonable cuts to discriminate signal from 
background and presenting the kinematics of the events. We also discuss background 
subtraction. We focus on the case of massive Xi particles and determine with what 
accuracy A might be measured. 

2.2 The Nature of the LSP 

Supersymmetric collider signals are notoriously complicated and depend strongly 
upon the parameter space, even if one restricts oneself to the MSSM. The monojet 
signal of interest for this paper depends primarily on two properties of the supersym- 
metric model: (i) the precise nature of the LSP (e.g., the relative contributions of 
the bino, wino and higgsino components of the neutralino LSP wave function); and 
(ii) the branching ratio of the squark into a quark and the LSP. 

In i?-parity-conserving supersymmetric models, the most likely candidate for the 
LSP (excluding the gravitino which is not relevant for this discussion) is the lightest 
neutralino Xi- 4 I* 1 general, Xi is a linear combination of bino, wino and higgsino 
interaction eigenstates. The relative contributions of each of these components de- 
pends on the parameters that govern the neutralino mass matrix. These include mz, 
the higgsino mass parameter /x, the gaugino mass parameters Mi, M 2 , and the ratio 
of Higgs vacuum expectation values tan /3. (Approximate formulae for the neutralino 
masses and mixing matrix in terms of these parameters can be found in Ref. [45].) 
If the higgsino component were dominant, then the monojet signal proposed in this 
paper would not be viable, as the corresponding cross-section for squark-neutralino 
production would be suppressed by a light quark Yukawa coupling. In this case, other 
methods must be employed to test for supersymmetric coupling relations. However, 
we note that in many (though not all) models of supersymmetry, the parameter 
is parametrically larger than Mi, M 2 and mz, in which case the higgsino component 
of the LSP is small. Henceforth, we will assume that the underlying SUSY model 
supports a Xi th & t is dominantly gaugino in nature. 

This still leaves the question of the relative contributions of the bino and wino 
components of the Xi wave function. In any supersymmetric model with the unifica- 
tion of tree-level gaugino mass parameters, the ratio of the low-energy values of the 

4 Cosmological and laboratory constraints rule out other possible candidates such as the sneu- 
trino [42], gluino [43] or charged slepton [44] in almost all possible parameter regimes. 
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gaugino mass parameters is given by [46], 

5a' 2 

Mx ~ -^-M 2 ~ 0.5M 2 . (2.1) 

Assuming that the tree-level gaugino masses are non- vanishing (and are of order mz), 
then the bino component of the Xi wave function is dominant. This is typical of most 
mSUGRA models [46], but is more general and depends only on the assumptions 
outlined above. 

Alternatively, it is possible that the tree-level gaugino masses vanish, in which 
case eq. fl2-l|) (which holds trivially) is irrelevant. In this case, the gaugino mass 
parameters arise at one-loop. In particular, a model-independent contribution to the 
gaugino mass is present whose origin can be traced to the super-conformal (super- 
Weyl) anomaly, which is common to all supergravity models [47]. This contribu- 
tion is dominant in models of anomaly- mediated supersymmetry breaking (AMSB). 
Eq. (|2.1|) is then replaced (in the one-loop approximation) by 

^ - , (2.2) 

where m 3 / 2 is the gravitino mass (typically assumed to be on the order of 1 TeV), 
and the bi are the coefficients of the MSSM gauge beta-functions corresponding to 
the corresponding U(l)y, SU(2)l, and SU(3) gauge groups: (bi,b 2 ,b 3 ) = (y, 1, —3). 
Eq. fl2.2p yields Mi ~ 2.8M 2 , which implies (under the assumption that is some- 
what larger than M 2 ) that the lightest chargino pair and neutralino comprise a nearly 
mass-degenerate triplet of winos over most of the MSSM parameter space. Typically, 
the corresponding neutralino is the LSP, whereas the wino is the next-to-lightest su- 
persymmetric particle (NLSP). 5 

The squarks that are produced in Fig. [I] are primarily first-generation squarks, 
either q^ or q^ (as the mixing between these states is negligible). In the case where 
the LSP is dominantly bino-like, then the branching ratio B(<£r — > qxi) — 100%. In 
contrast, the dominant channel for qi decay is into the heavier wino-like neutralino 
which subsequently decays into the LSP. Thus, the production of qn in association 
with the LSP is more likely to produce a monojet. In the case where the LSP is 
dominantly wino-like, the reverse is true, and the production of q^ in association 
with the LSP is more likely to produce a monojet. 

In this paper we shall examine several scenarios. The first scenario considered 
will be an optimistic mSUGRA point, which will illustrate the effect of the SUSY 
backgrounds, making a coupling extraction from LHC data difficult. Next, we shall 
assume a wino dominated neutralino, where signal cross-sections are higher and 

5 We ignore the gravitino, which could be lighter than the neutralino. Nevertheless, in almost all 
cases, the neutralino behaves as if it is the LSP. That is, the neutralino will be sufficiently long-lived 
so that it will always escape the collider detector before it decays. 
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easier to measure over the background. In this case, the coupling extraction from 
LHC data is somewhat easier. If the lightest neutralino is light and found to be wino 
dominated, then it does not constitute a significant portion of a thermal dark matter 
relic, since winos annihilate too efficiently in the early universe [48]. 

3. Testing the Supersymmetric Coupling Relations with Mono- 
jets 

3.1 The Signal Process 

We define our parton-level signal processes in the following way: 

• For a bino LSP: g + u — > Xi + ur, g + d — > Xi + d>R as displayed in Fig. [l] 
followed by u R — >■ uxl or <i R — > c?%°. Here, A = g r , the U(l)y gauge coupling 
of the SM. The charge-conjugated processes are also included. 

• For a wino LSP: g + u — > Xi + ul, g + d — > Xi + d>L, 9 + u Xi + di, 
9 + d ->■ Xi + u L , where u L uxi/dxt, d L dxi/uxi and xl -» S + x?, 
where S is either QCD radiation too soft to be identified as a jet, or a lepton so 
soft that it passes the lepton veto. For wino LSPs, the mass splitting between 
Xi and x? is ~ 200 MeV, so the dominant decay xX ~ Xi n+ typically results 
in a soft pion that is too soft to be measured by the LHC experiments. Both 
the Xi an d xt are wino-like and are related via a SU(2)l transformation, so the 
production process amplitudes are constrained by the MSSM to be proportional 
to A, equal to the SU(2) L gauge coupling g. The charge-conjugated processes 
are also included. 

In principle, processes involving heavier non-valence quarks like g + c — > Xi + c 
also contribute to our signal and SUSY backgrounds. We have only included initial 
states with a g,u,u,d or d, since processes with quarks beyond the first generation 
are negligible. Since we assume that Xi is stable, it leaves no direct trace in the 
detector, and so the signal consists of a jet recoiling against apparently missing 
transverse momentum of magnitude jj T . In general, j/ T is measured to be different 
to the transverse momentum of the jet pr(i) because of measurement errors and also 
because of soft QCD radiation, which may not be included in the jet, but measured 
in the calorimeter nonetheless. 

3.2 Major Backgrounds and Basic Cuts 

As we have seen in the last section, our signal consists of a hard jet recoiling against 
missing energy. Therefore, the major SM backgrounds are [26,49] 

• Z(— > v z/)+jet. Z production in association with a hard jet, where the Z decays 
into a pair of neutrinos. 
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• W(— > Tis)+jet. W plus jet production, where the W decays into a tau and a 
neutrino, and the tau is either not detected, or lost in the jet. 

• W(— >■ e/pL z/)+jet. W plus jet production, where the W decays into an electron 
or a muon and a neutrino and the electron/muon is undetected or lost inside 
the jet. 

• QCD jet production together with mismeasurement of the energy deposited in 
the detector. One could produce di-jets, for instance, and one of the jets could 
be lost in the detector (or its energy mismeasured so that it fluctuate below 
the transverse momentum required to identify the jet). 

In principle, di-vector boson production VV , with V, V = W, Z is an additional 
source for the SM backgrounds if one vector decays into a neutrino and the other into 
jets. However, the cross section is much smaller than Z/W+]et and VV production 
can be safely neglected in our analysis [28]. 

In order to reject most of the W+jet background, we employ a lepton veto. We 
veto events with an isolated electron or muon with p? > 5 GeV and with |?7| < 2.5. 
The isolation criterion demands less than 10 GeV of additional energy in a cone of 
radius AR = a/ A0 2 + Arf = 0.2 around the lepton momentum, where and r] are 
the azimuthal angle and pseudorapidity, respectively. 

After the lepton veto, Z(— > uu)+]et remains the most important irreducible 
background. Fortunately, this background can be directly derived from data itself 
by measuring Z(— > e + e"/ / u + /i~)+jet and relating both processes with the help of the 
measured Z branching ratios. Systematic uncertainties are in this case significantly 
reduced [50]. However, this comes with the cost of higher statistical uncertainties, 
because the Z(— > z/z/)+jet cross section is roughly three times larger than the Z[— > 
e + e~ / fx + /i~)+jet cross section. 7+jet cross-sections at high p T could be used in the 
future to estimate the Z(— > uu)+]et cross section with better statistics. 

We will use two statistical estimators; cf. Refs. [26,49]. The first optimistic esti- 
mator takes only the statistical fluctuations from Z{— >■ uu)+]et into account. In this 
case the significance is given by S/ VB with S (B) the number of expected signal (SM 
background) events. This case would apply when the background Monte-Carlo is so 
well tuned and tested with LHC data that its output may be fully trusted. We also 
employ a conservative estimator, where the statistical fluctuations are dominated by 
the Z(— > e + e~ / /j,~)+jet calibration sample. According to Refs. [26,49], this case 
corresponds to a significance of S/ yJB. & The conservative estimate would be used 
in the case that one does not trust at all the Monte-Carlo background calculation, 
and instead measures the Z+jet background from LHC data. We expect the true 

6 In Refs. [26, 49] the transverse momentum of the monojct was required to lie above 1 TeV. 
In contrast, we will use less hard cuts for which the Z(— >• e+e~/,u + ^~)+jct acceptance is slightly 
better; see Ref. [49]. However, we will always use S/VYB as a conservative estimator. 
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(i.e. measurable) significance to lie between our two estimators after the Monte Carlo 
estimates [50] have been properly tuned to LHC data, including W+jet/7+jet. The 
modeling of systematic error on the measurements of SM backgrounds is an experi- 
mental question, and as such is beyond the scope of this paper. 

After cuts, once the large Z(— > vP)+jet background is subtracted, we will still 
need to subtract the VF+jet background, which is larger than our signal. This can be 
measured in the case where the lepton from the W is visible, by extrapolating into 
the region where the lepton is invisible, either because it is lost in the jet (in which 
case one could extrapolate in AR between jet and lepton), or because the lepton is 
missed. 

An estimate of the QCD background can be found by full detector simulation [50]. 
The quantity of interest for estimating these backgrounds is the jet energy response 
function (JERF) R, equal to the ratio of the measured jet energy to the true one. 
We fit the JERF to the full detector simulation results in [50] in the same spirit 
as Ref. [51]. In our analysis, R is used to scale all of the components of the jet 
four-momentum. R is well fit by a probability distribution function 

p(R) = -^2L + 0.01 A e 7 - 32R 0(0.9 - R), (3.1) 

V27T(7 2 

where is the Heaviside step function, A is an unimportant normalization constant 
(« 0.01), and 

" 0.03 (3.2) 



y/E/GeV 

is the resolution of the Gaussian part of the JERF, and © denotes the fact that 
the two terms are added in quadrature. The second term in p(R) reflects the 1% 
probability that a significant portion of the jet will go unmeasured due to cracks 
in the detector and other effects. This exponential tail is responsible for the QCD 
background to our SUSY monojet signal, and in order to increase statistics on the 
sample, we impose that either the hardest or second hardest jet in the sample be in 
the Gaussian tail, taking into account the additional factor of 100 needed to calculate 
the cross-section after cuts. Other jets are all drawn from p(R). QCD backgrounds 
to monojets will be extrapolated from data, for example from di-jet, 7+jet and 
"Mercedes" type 3-jet events [50], although we note that these techniques may not 
address more conspiratorial backgrounds, which could require the use of tracking 
information. QCD backgrounds will then be subtracted from event samples, leaving 
statistical fluctuations only once the systematics have been dealt with. The JERF 
was only applied to QCD backgrounds, not to the other samples, since it should 
only have a small effect upon them. We note that Herwig++ includes b backgrounds 
within this QCD sample, so the case of the production of bb, where one of the bottom 
quarks' momenta primarily goes into a neutrino resulting in a monojet signature, is 
included in our estimate. 



- 11 - 



A further background source is supersymmetry itself. For example, squark pair 
production and the subsequent decay of the squarks into the Xi LSP and a quark is a 
possible SUSY background. If the two jets from the squark decays overlap they can 
appear as a single monojet. Initial-state radiation on X1X1 production also provides 
a SUSY background, although this is usually rather small. Note that the background 
of squark pair production followed by the decay into two jets plus two X?, where one 
of the jets is in the tail of the JERF, is not in our SUSY background sample. We 
expect this background to be smaller than the QCD background, which is already 
small. 

In the case of a bino-like LSP, X2 + Ql production followed by xl ~ * Xi + v ^ and 
qi, — > x? + q would also produce monojet events. The cross-section for X2 + Ql is 
enhanced by Kg 2 / g' 2 ~ 2 over xi + Qr production (where the kinematical factor K 
provides a suppression due to the heavier x° mass). However, due to the branching 
ratio factors for the final state decays, the number of x° + Ql events that are observed 
as monojets is significantly smaller than our signal and can be neglected. 

In the case of a (not too heavy; see Section |5.4j) wino-like LSP, the dominant 
SUSY background is wino pair production plus a jet from initial state radiation. The 
cross section is quite large because two winos, WW, can be produced via a Drell-Yan 
like process [52,53], i.e. PP ->■ 7*/^* -> WW. Here, 7* (Z*) denotes a virtual 7 
(Z). In the wino-LSP scenario, both W means either Xi or Xi> both of which are 
approximately winos. Because the Xi is quasi-degenerate in mass with the x?> ft 
decays into a soft pion and the x^. The pion is typically too soft to be detected [48] 
unless special analysis techniques are used [54] and so the Xi is effectively invisible 
to the detector. The additional jet is then produced via initial state radiation. 

We have employed Herwig++2.4.2 [55-57] 7 to simulate the signal and the back- 
grounds at tree level in pp collisions at a center-of-mass energy of 14 TeV. An overview 
of the simulated signal and background samples is given in Table [TJ. We have in- 
cluded the pair production and the two and three-body decay of all SUSY particles. 
All events from this sample that yield the correct monojet topology that are not 
our signal process are classified as SUSY background. The Z+jet/IU+jet back- 
grounds are also obtained with Herwig++. The Herwig++ output was analyzed using 
HepMC-2.04.02 [59] and ROOT [60,61]. Jets were reconstructed with the help of 
f astjet-2.4. 1 [62]. We employed the anti-fct jet algorithm with R = 0.7 [63]. Only 
jets with p T > 30 GeV are used in our analysis in order to be less sensitive to un- 
derlying event modeling. In Table [j], we list the total number of simulated events 
for each sample. In the bottom half of the table, the total number of SUSY pair 
production events is listed first, with the number in the subset corresponding to our 
signal events listed under "comments" . 

7 We have used a modified version of Herwig++2.4.2 which is also able to deal with negative 
(Majorana) gluino masses; see the following link for details [58]. 
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sample 


simulated events 


comments 


Z[— > uu)+]et 
W(-t ez/)+jet 
W(-> fiv)+iet 
W(-t rz/)+jet 
QCD 


2140 000 
2 960 000 
2 960 000 
2 960 000 
18 650 000 


pr(jet) > 70 GeV 
Pr(jet) > 70 GeV 
Pr(jet) > 70 GeV 
p T (jet) > 70 GeV 
Pr(jet) > 30 GeV 


mSUGRA 
mAMSB 
others 


3 280 000 
2 940 000 
100 000 


Included 10 300 signal events. 
Included 46 700 signal events. 
Number of signal events varies. 



Table 1: Monte Carlo samples of SM backgrounds and SUSY events used in our analysis. 
All simulated events were generated with Herwig++2 .4.2 for pp collisions at a center- 
of-mass energy of 14 TeV. The Z+jet and VU+jet samples correspond to an integrated 
luminosity of 5 fb _1 . The QCD sample corresponds to 30 nb _1 , the mSUGRA sample 
corresponds to 20 fb _1 and the mAMSB sample corresponds to 100 fb _1 . The cut on 
pr(jet) is performed at the parton level. The mSUGRA [mAMSB] sample corresponds 
to sparticle pair production assuming the benchmark scenario eq. (|4.1|) [eq. (fTj])]; the 
number of SUSY events that correspond to our signal monojet events is indicated in the 
third column of the Table for these two rows. The scenarios denoted by others are those 



of the parameter scan of Section 5.4; the corresponding integrated luminosities have been 



adjusted for each scan point such that 100 000 events for each point are generated. 

4. Measurement of the }£{qq Coupling for a Bino LSP 

In this section we will show that it is possible to estimate the XiQrQ coupling for the 
case of a bino-like Xi LSP from data. We choose an mSUGRA benchmark scenario 
for investigation. 



4.1 Benchmark Scenario 

Our benchmark scenario is a light mSUGRA scenario 
It is described by the parameters 



with a binodike Xi LSP. 



M = 220 GeV, M 1/2 = 180 GeV, A = -500 GeV, tan = 20, sgn(jx) = +1 . (4.1) 

For this scenario, the XiQrQ coupling is given by A = 0.99g', where the deviation from 
A = g' is due to the small admixture of wino and higgsino eigenstates in the LSP. All 
spectra are calculated with S0FTSUSY3 . . 13 [65] and fed into the event generator 
via the SUSY Les Houches Accord [66]. The sparticle masses are given in Table ||. 

This scenario has a relatively light SUSY mass spectrum and a total signal cross- 
section of a(pp qnXi) = 520 fb. For spectra a little heavier, the monojet cross 
section is too small and could not be seen above the SM backgrounds. Note that the 
mSUGRA scenario, eq. ( |4.1| ), lies at the edge of the region excluded by the Tevatron, 
but is still allowed [67,68]. 
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sparticle 


mass [GeV] 


sparticle 


mass [GeV] 


xS 


70.2 


xl 


365 


xi 


132 


xi 


370 


xl 


133 


h 


378 




189 


b 2 


443 


tl 


226 


ur/c r 


454 




230 


d R /s R 


455 




234 


9 


456 




242 


u L /c L 


463 


eL/J^L 


255 


d L /s L 


470 


r 2 


259 


t 2 


477 


~o 

X3 


359 







Table 2: Sparticle mass spectrum of the mSUGRA benchmark scenario given by eq. (|4.1|). 
The sparticles are ordered by their mass. 

4.2 Event Numbers and Cuts 

In this section, we will develop a set of cuts that allow a measurement of the XiQrQ 
coupling assuming our mSUGRA benchmark scenario, eq. ( |4.1|) . On one hand, the 
number of signal events needs to significantly exceed the statistical fluctuations of the 
SM backgrounds. On the other hand, we desire a good signal to SUSY background 
ratio if we want to measure A to a high precision. 

In Table |3|, we display the total number of SM background events (second col- 



cut 


all SM SUSY bkg. signal 


s/Vb 


trigger 
lepton veto 
number(jets)=l 
tf T > 180 GeV 
m(jet) < 70 GeV 


1.14 x 10 8 2.91 x 10 7 130 000 
7.57 x 10 7 1.76 x 10 7 130 000 
3.35 x 10 7 55 900 35100 
3.28 x 10 6 32 300 22 300 
3.00 x 10 6 12100 20100 


6.1 (2.3) 
12 (4.7) 
12 (4.4) 


tau veto 
6-jet veto 


2.75 x 10 6 9 950 20 000 
2.66 x 10 6 9 290 20 000 


12 (4.6) 
12 (4.6) 



Table 3: Cut flow for the mSUGRA benchmark scenario in eq. ( fl.l|) . We present 
the cuts in the first column and the number of SM, SUSY background and signal events 
in the second, third and fourth column, respectively. We also show in the fifth column 
the resulting significance for the monojet signal, i.e. B corresponds to the number of SM 
background events and S is the number of signal events. The significances in brackets 
are our conservative estimate, i.e. S/VYB. We have assumed an integrated luminosity of 
300 fb _1 at -y/i = 14 TeV. Note that in addition to the j/ T cut, the same cut has been 
applied on the jet-pr- 
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umn), the number of SUSY background events (third column) and the number of 
signal events (fourth column) for different sets of cuts (first column). We also show in 
the fifth column the resulting significance for the monojet signal, i.e. B corresponds 
to the number of SM background events and S is the number of signal events. Note 
that B refers to the SM background only, but the SM contribution completely domi- 
nates the background. We assumed an integrated luminosity of 300 ftT 1 at y/s = 14 
TeV. We cannot give meaningful numbers for the significances in the first two rows 
because SM backgrounds additional to those of Section ^]2| would need to be included. 

As a first cut (denoted by trigger), we demand at least one jet with transverse mo- 
mentum, pr(jet), larger than 100 GeV. In addition, the amount of missing transverse 
momentum, j/ T , must also exceed 100 GeV. These cuts correspond to the planned 
one-jet plus missing energy trigger used by the ATLAS collaboration [69] for the 14 
TeV run. Therefore, all events in the first row of Table |] will be recorded. Note 
that CMS plans to use harder cuts [70]. We can see in Table |] that two orders of 
magnitude more SUSY background events will pass the trigger than signal events. 

The SUSY background is reduced by roughly a factor of two after we apply (in 
addition to the trigger cut) a veto on isolated electrons and muons; see Section |3~2| . At 
the same time, the number of signal events is nearly unaffected. The SM backgrounds 
are also reduced by the lepton veto, because most of the W+jet background events 
will not pass this cut. However, the Z(—> up)+]et background stays nearly the same. 

To reduce the SUSY background further, we also apply a veto on a second 
jet if its pt is larger than 30 GeV and \rj\ < 5. This corresponds to a veto on a 



second jet, because we only count jets above 30 GeV, cf. Section |3]2|. The number of 
SUSY background events now has the same order of magnitude as the signal, namely 
(9(10 4 ) events. Although the signal possesses no second jet at parton level, we might 
produce one due to initial and final state radiation, as is borne out by the third row 
of Table [| In particular, a veto on a second jet also reduces the number of signal 
events by roughly a factor of four. 

At this stage we are able to give some reliable numbers for signal significances, 
assuming that backgrounds can be reliably constrained and subtracted. The numbers 
in brackets correspond to our conservative estimate as described in Section i.e. we 
assume that the dominant Z(—> uu)+jet background is estimated purely from a 
measurement of Z(—> e + e~//! + /i - )-|-jet. We observe an optimistic (conservative) 
significance for the monojet signal of 6.1 (2.3). 

In order to increase the signal to SM background ratio we make use of the jf T 
and jet-pr distribution of our monojet signal. They possess a Jacobian peak, where 
the position depends on the qn and Xi masses (see Appendix A). We present the jf T 
distributions in Fig. || The Jacobian peak of our monojet signal (red histogram) 
lies around 180 GeV. At the same time the distributions of the SM backgrounds 
(black, blue, magenta and dashed black histograms) fall off exponentially. An ad- 
ditional lower bound on jf T and the monojet of 180 GeV, raises the optimistic 
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pTmiss [GeV] 

Figure 3: Missing transverse momentum distribution for the Z(—> vD)+jet (black his- 
togram), W{— > r^)+jet (blue histogram), W{— > e^)+jet (magenta histogram), the SUSY 
background (green histogram), the QCD background (black dashed histogram) and the 
signal process (red histogram). The first three cuts of Table || have been applied. The 
number of events correspond to an integrated luminosity of 100 fb _1 at -y/s = 14 TeV. We 
assume the mSUGRA scenario of eq. ( |4.1| ). 

(conservative) signal versus SM background significance to 12 (4.7). The required 
level of cancellation is high: around 1 in 100 for a high jj T cut. This means that the 
experimental systematic on the dominant background also needs to be at the level of 
< 1%, a challenging (but not obviously impossible) proposition. The experimental 
systematic on the subdominant PU+jet background subtraction would need to be at 
the level of a few percent. 

Our goal is to reconstruct the XiQrQ coupling as precisely as possible. In order 
to accomplish this goal, we need to reduce the SUSY background further while 
leaving the number of signal events unchanged (in order to have a visible signal). 
For this purpose, we have examined the invariant mass of the hardest jet, m(jet), 
which is presented in Fig. |] for the signal (red histogram), the SUSY background 
(green histogram) and the Z+jet background (black histogram). Each histogram is 
separately normalized so that the total number of events is one. The first four cuts of 
Table [3] are applied. Note that the shape of m(jet) for VU+jet and QCD backgrounds 
follows those of Z+jet. 

We first observe that the distribution of the Z+jet background looks very similar 
to the signal distribution. This is expected; the m(jet) distribution for these two 
cases is what is predicted by an approximately massless initial parton (which ideally 
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140 160 180 
m(jet) [GeV] 

Figure 4: Invariant mass distribution of the hardest jet for the Z(—> uP)-\-jet background 
(black histogram), the SUSY background (green histogram) and the signal process (red his- 
togram). The first four cuts of Table || have been applied. The distributions are normalized 
to one. For the signal, we assume the mSUGRA scenario of eq. (|4.l| ). The TU+jet and 
QCD backgrounds have a distribution almost indistinguishable from that of Z(—> vi>)+jet. 



would peak at zero mass) when one applies lower py(jet) cuts, cutting out the very 
low mass region. Therefore, a cut on m(jet) cannot increase the signal significance 
of the signal over SM background. However, the SUSY background distribution 
has a very different shape. It is relatively flat compared to the signal and has its 
maximum at a larger m(jet) value. The (on average) larger jet invariant mass of the 
SUSY background stems mainly from events, where a pair of squarks is produced and 
where the jets from squark decays go roughly in the same direction and are identified 
as only one jet. Relative to the cuts listed in Table |3], imposing the requirement that 
m(jet) is less than 70 GeV reduces the SUSY background by roughly a factor of 
three and leaves the number of signal events nearly unchanged. 8 Note that the 
cut on m(jet) mostly suppresses the high-j/ T SUSY background events in Fig. [| 
i- e - tIt ~ 300 GeV. In the next section, we will employ these cuts, i.e. the first five 
cuts in Table ^|, for the coupling reconstruction. 

Finally, in the last two rows of Table [3] we show the effects of a tau-lepton 
veto and a subsequent 6-jet veto. In both cases we assumed an identification (ID) 
efficiency of 100%. We regard a tau as identified if it decays hadronically and has 
Pt > 15 GeV and \r]\ < 2.5. We observe that neither veto helps much, and so their 
effect would be even more reduced if realistic ID efficiencies were assumed. We see 

8 For the mSUGRA benchmark point, the number of surviving monojet events after the cut on 
m(jet) arising from x§ + production (which is included in the SUSY background in Fig. is 
about 5% of our signal events. 
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that the SUSY background would be reduced by 18%. However, the ID efficiencies 
are doubtless too optimistic and in reality the suppression of the SUSY background 
would lie below 11% [50,71]. Thus, we will not employ the tau and 6-jet vetoes here. 

4.3 Coupling Reconstruction 

We now show to which precision the XiQrQ coupling A of the mSUGRA benchmark 
scenario of Section [4.1| can be reconstructed with an integrated luminosity of 300 fb _1 
at y/s = 14 TeV. We employ the cuts developed in Section [O] (without the tau and 
6-jet veto). Since the signal monojet cross section is proportional to A 2 , the relative 
error on A will be roughly one half of the relative error on the monojet production 
cross section. Note that we reconstruct A under the assumption that the SUSY 
backgrounds come from our mSUGRA point. 



error 


A^mono /'-'"mono 


AA/A 


luminosity 


3.0% 


1.5% 


PDF uncertainty 


17% 


8.3% 


NLO corrections 


16% 


8.0% 


sparticle mass Am = 10 GeV 


11% 


5.6% 


statistics (optimistic) 


8.6% 


4.3% 


statistics (conservative) 


23% 


11% 


total (optimistic) 


27% 


14% 


total (conservative) 


35% 


17% 



Table 4: Relative errors for the signal monojet cross section (second column) and the 
XiQrQ coupling A (third column) from different sources (first column). The numbers are 
for the mSUGRA benchmark scenario, eq. ( [4.11 ). 

In Table f|, we show the statistical and the most important systematic errors for 
a measurement of the signal cross section and XiQrQ coupling A. The total error is 
obtained by combining all errors in quadrature. The errors were estimated as follows: 

• In order to determine the monojet cross section from the observed number of 
events, one needs to know the total number of collisions, i.e. the integrated 
luminosity. According to Ref. [72-74], the integrated luminosity at LHC will 
be measured to a precision of 3% or better. 

• To translate the hadronic cross section to a parton-level cross section, one needs 
to know the parton distribution functions (PDFs) accurately. We estimated 
the error from the PDFs by comparing the hadronic cross sections with differ- 
ent PDF releases; namely CTEQ61, CTEQ611, CTEQ6m [75], MSTW20081o, 
MSTW2008nlo, MSTW2008nnlo [76], GJR08 [77,78] and a cubic interpolation 
of MRST LO** [79] (default Herwig++ PDF). We found the largest variation 
between MRST LO** and MSTW20081o. It corresponds to a variation of ±17% 
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around the central value of the total hadronic signal cross section. We expect 
the PDF errors to be reduced significantly after the input to PDFs of LHC 
data, but for now we use the current value of the error. 

• Further uncertainties arise from the unknown next-to-leading order (NLO) cor- 
rections, which are not yet fully implemented [80]. The corresponding uncer- 
tainty is estimated by varying the renormalization scale Q between \{rriq R + 
m^o) < Q < rriq R + m^o (where the average squark mass of the right-handed 
first and second generation squarks is denoted by Triq R ) leading to an error of 
+16%/ — 13% for the total LO cross section. In Table 0, we use 16% in order to 
be rather conservative, with the knowledge that the NLO calculation would be 
ready and significantly decrease (possibly by a factor of a few) this uncertainty. 
The LO cross section was derived with Prospino2. 1 [81]. 

• The relatively light SUSY spectrum in Table || leads to copious production of 
SUSY particles at the LHC [82]. In this case, we expect that the x? and squark 
masses can be reconstructed to a precision of at least 10 GeV [4]. Varying the 
squark and Xi masses by ±10 GeV introduces an error of ±11% around the 
central value of the total monojet cross section. 

• We always get an error due to statistical fluctuations of the SM backgrounds. 
In Table 0, we present two statistical estimates for this error as described in 
Section |3.2| . The optimistic (conservative) estimator is \JB (v7 B), with B the 
number of SM background events. 

We have neglected in Table |] errors from the SUSY background, because an pre- 
cise estimate of its systematic uncertainties lies beyond the scope of this publication. 
However, we expect them to be small, because it is in principle possible to extrap- 
olate them from data. As we observe in Fig. |], the invariant mass distribution of 
the monojet is dominated by the SUSY background for m(jet) > 70 GeV providing 
a control sample for a fit of the SUSY background. Further model assumptions will 
also reduce the systematic uncertainties. 

We can see in Table ^ that the biggest uncertainties are: the PDF uncertainty, 
NLO corrections and from the SM statistics (conservative estimate). However, one 
can hope that the uncertainty of the PDFs will rapidly decrease after LHC has taken 
and analyzed some first data. Furthermore, the calculations of higher order correc- 
tions might improve in the future. Unfortunately, an integrated luminosity of 300 fb _1 
corresponds already to the maximum expected LHC data. Some improvement might 
be possible by combining the data of both multi purpose detectors. 

To conclude, a measurement of A via monojet production at the LHC is in prin- 
ciple possible for scenarios with a bino-like Xi LSP like mSUGRA. However, at least 
parts of the SUSY mass spectrum need to be relatively light, i.e. the right-handed 
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up-squark and down-squark (xi LSP) mass should be < 500 GeV (< 100 GeV). A 
measurement of the XiQrQ coupling at the 10%-level is feasible by the end of LHC 
running. The bino couplings to quark-squark pairs is proportional to the correspond- 
ing squark hypercharge Y, which is different for cIr, and ur. Since production of both 
squarks contributes to our SUSY monojet signal, this fact must be taken into account 
when extracting a value of A from the data. 9 For much larger masses, the monojet 
cross section is too small to allow its extraction. 

The situation is much more promising for SUSY scenarios with a wino-like Xi 
LSP. The monojet cross section is now enhanced compared to the bino-LSP case due 
to larger gauge couplings and due to additional diagrams involving a charged wino. 

5. Reconstruction of the )$qq Coupling for a Wino LSP 

5.1 Benchmark Scenario 

The second benchmark scenario is a mAMSB scenario [47,48,83] with a wino-like Xi 
LSP and with a wino-like xt NLSP that is nearly degenerate in mass with the Xi 
LSP. The scenario is described by the parameters 

M 3/2 = 33TeV, M = 200 GeV, tan/3 = 10, sgn(/i) = +1 . (5.1) 

The resulting sparticle masses at the electroweak scale are given in Table |5]. 



sparticle 


mass [GeV] 


sparticle 


mass [GeV] 


X? 


106.5 


xl 


593 


xi 


106.7 


xi 


594 


Tl 


113 




634 


V T 


135 


h 


688 




138 


u L /c L 


722 




150 


J 2 


723 




159 


d L /s L 


726 


r 2 


179 


ur/c r 


726 


~o 
Xl 


298 


d R /s R 


732 




521 


9 


745 


~o 
Xi 


584 







Table 5: Supersymmetric particle masses for the mAMSB scenario defined in eq. (|5.l|). 
We show the masses of the wino-like Xi an d xi U P to a precision of four digits, because 
they are nearly degenerate in mass. The sparticles are ordered by mass. 

9 We expect that the first generation squarks are roughly mass-degenerate. However, if there were 
a significant hierarchy between the up and down-type squark masses, then one would dominantly 
produce the lightest squark, and the extraction of A from data would be more straightforward. 
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In this scenario, the total signal cross-section is a{pp — > qiXi/xt) — 470 fb, and 
A = 0.99^ due to the small bino/higgsino admixture. The lightest chargino, which 
will be included in our signal, couples with strength 0.98(7 to d L u because of a small 
higgsino admixture and with l.OOg to u^d. The small differences between the xt an d 
Xi couplings are much smaller than our measurement errors, and shall be neglected 
in the following. 

5.2 Event Numbers and Cuts 



In this section we apply the same procedure as we did in Section [Q[ We develop a set 
of cuts such that a clear signal over the SM backgrounds is visible. Furthermore, we 
want to obtain a good signal to SUSY background ratio, which is needed for a precise 
estimation of the XiQlQ coupling A. There are some differences in our analysis to 



the bino-LSP scenario, eq. (|4.1| ). Instead of an integrated luminosity of 300 fb , we 
assume only 100 fb -1 to show that less statistics is needed to reconstruct A. This is 



due to the much larger monojet cross section for wino-LSP scenarios; cf. Section [IT 
For the same reason the discovery potential is better, i.e. we can investigate a heavier 
spectrum; see Table || and Table |. 

The cut flow for the mAMSB benchmark scenario, eq. ( |5.1|) , is given in Table |[ 
The first two cuts are identical to those in Table [| However, the jet veto is slightly 
different (third row in Table ^). Instead of vetoing additional jets with pxs larger 
than 30 GeV, we now relax this cut to 50 GeV. 



cut 


all SM SUSY bkg. signal 


s/Vb 


trigger 
lepton veto 
p T (jet2) < 50 GeV 
rf T > 300 GeV 
m(jetl) < 80 GeV 


3.81 x 10 7 1.04 x 10 6 44100 
2.52 x 10 7 621000 43 800 
1.73 x 10 7 111000 16 200 
171000 11000 8 390 
135 000 6 020 6 370 


3.9 (1.5) 
20 (7.7) 
17 (6.5) 


tau veto 
6-jet veto 


119 000 5 840 6 370 
115 000 5 290 6 320 


18 (7.0) 

19 (7.0) 



Table 6: Same as Table fj, but now for the mAMSB benchmark scenario, eq. (5J-). We 
have assumed an integrated luminosity of only 100 fb -1 at -y/s = 14 TeV. Jetl and jet2 
denote the jet with the largest pr and second-largest px, respectively. 



In contrast to the mSUGRA scenario, the SUSY background for our mAMSB 



benchmark point is dominated by Drell-Yan wino pair production; cf. Section ^2 
At the parton level, the event topology is similar to that of the signal and to that 
of the dominant SM background: two invisible particles recoiling against a single 
quark or gluon. After initial and final state radiation, additional jets can arise and 
the event might fail the jet veto. The veto then reduces the relative number of signal 
and (Z+jet and wino-pair plus jet) background events equally due to the similar 
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Figure 5: Missing transverse momentum distribution for the Z(—> vD)+jet (black his- 
togram), W{— > r^)+jet (blue histogram), W{— > e^)+jet (magenta histogram), the SUSY 
background (green histogram), the QCD background (black dashed histogram) and the 
signal process (red histogram). The first three cuts of Table ^ have been applied. The 
number of events correspond to an integrated luminosity of 100 fb _1 at -y/s = 14 TeV. We 
assume the mAMSB scenario of eq. ( |5.1| ) . 



event topologies. Therefore, the jet veto reduces the significance S/ yB whereas the 
signal to SUSY background ratio stays roughly constant. A weaker jet veto thus 
increases the visibility of the signal. Note that for the mSUGRA scenario, eq. (|4.1|) , 
the topology of most of the SUSY background is different. For example, in the 
mSUGRA scenario, the dominant background is squark pair production, where both 
squarks decay into a quark and a x*?. 

After the jet veto, we still observe in Table |6| an overwhelming SUSY background, 
i.e. the signal to SUSY background ratio is roughly seven whereas for our mSUGRA 
benchmark point this ratio was less than two; see Table [5]. However a great improve- 
ment is possible with a cut on the jj T . The relevant distributions (after the first 
three cuts in Table |6| were applied) are shown in Fig. |5|. Here, the SUSY background 
distributions as well as the SM background distributions fall off exponentially. In 
contrast, the j/ T distribution of the SUSY background in Fig. |] does not show the 
same exponential fall-off; a clear sign that the SUSY background is different in the 
mSUGRA and mAMSB cases. Applying again a lower pj T cut around the Jacobian 
peak of the signal (the red histogram in Fig. [5]) increases the signal to SUSY back- 
ground ratio to 0.8 and the optimistic (conservative) significance to 20 (7.7). Note 
that we have also applied the same cut on the px of the hardest jet. We also see 
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Figure 6: Invariant mass distribution of the hardest jet for the the Z+jet background 
(black histogram), the SUSY background (green histogram) and the signal process (red 
histogram). The first four cuts of Table ^ have been applied. The distributions are nor- 
malized to one. For the signal, we assume the mAMSB scenario of eq. ( |5.1| ). The QCD 
background has a distribution almost indistinguishable from that of Z{— > uP)+jei. 



that for this point, the required level of cancellation of the SM model backgrounds 
is at the ~ 10% (or larger) level. Deriving a systematic error smaller than this on 
the backgrounds measurement ought to be easily achievable with the high number 
of data expected. 

In Fig. H we present the invariant mass distribution of the hardest jet after the 
first four cuts of Table |] were applied. Because the SUSY background is dominated 
by Drell-Yan like wino pair production, the m(jetl) distribution is less flat than those 
in Fig. [|. However, an improvement with an upper cut on m(jetl) is still possible. 
Demanding m(jetl) < 80 GeV increases the signal to SUSY background ratio to 1.1. 
At the same time, the signal significance is not significantly changed. It is 17 (6.5) 
for our optimistic (conservative) estimate. We finally show in the last two columns 
of Table |]the effect of a tau and 6-jet veto. Again, we observe no great improvement 
by such a veto and thus we shall not apply it in the following. 

5.3 Coupling Reconstruction 

We now show to which precision the XiQlQ coupling A of the mAMSB benchmark 
scenario of Section |5.1| can be reconstructed with an integrated luminosity of 100 fb~ 
at y/s = 14 TeV. We employ the cuts given in Section ^]2] without the tau and 6-jet 



veto. We will closely follow the procedure employed in Section 4.3. We reconstruct A 



under the assumption of the SUSY backgrounds of our mAMSB point. We exhibit in 
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error 


^^"niorio / ^"mono 


AA/A 


luminosity 


3% 


1.5% 


PDF uncertainty 


17% 


8.3% 


NLO corrections 


18% 


9% 


sparticle mass Am = 10 GeV 


7.3% 


3.7% 


statistics (optimistic) 


5.8% 


2.9% 


statistics (conservative) 


15% 


7.7% 


total (optimistic) 


26% 


13% 


total (conservative) 


30% 


15% 



Table 7: Relative errors for the signal monojet cross section (second column) and the 
XiQLQ coupling (third column) from different sources (first column). The numbers are for 



the mAMSB benchmark scenario, eq. (5.1 



Table [I] the statistical and the most important systematic errors for a measurement 
of the cross section and thus the XiQlQ coupling. 

Note that errors from the SUSY background are not included. For our mAMSB 
benchmark scenario (and for many wino LSP scenarios in general), the SUSY back- 
ground is dominated by wino pair production plus jet. The differential distribution as 
a function of the missing transverse momentum is known [35]. One can also measure 
the SUSY background below the Jacobian peak of the monojet distribution. For our 
mAMSB benchmark scenario, we found that the SUSY background in an jf T interval 
of 150 GeV to 300 GeV is more than six times larger than the monojet signal as can 
be seen in Fig. |5|. At the same time, for the total number of SUSY events we found a 
(conservative) significance over SM background of roughly ten. Therefore, it should 
be possible to estimate the SUSY background from data which reduces significantly 
the systematic uncertainties of the SUSY background. Assumptions about the un- 
derlying model can further reduce this error. We thus expect that errors from SUSY 
background will not significantly alter the results in Table [7]. However, a precise 
reconstruction of the SUSY background lies beyond the scope of this publication. 



As in Section [O] the largest errors come from the PDF uncertainty, from the 
NLO corrections and from the SM statistics (conservative estimate). There is a good 
chance that the PDF uncertainty will decrease after LHC has analyzed some first 
data. Also the calculations of some of the unknown higher order corrections might be 
performed in the future. Finally, increasing the integrated luminosity would decrease 
the statistical error. In contrast to the bino LSP scenario, where we already assumed 
an integrated luminosity of 300 fb" 1 , the coupling reconstruction is possible earlier 
for our wino LSP scenario, because we only assumed 100 fb _1 of data. For example, 
assuming an integrated luminosity of 300 fb~ x and an error from the PDFs and from 
higher order corrections of 10%, the relative error of A would decrease from 13% (15%) 
in Table |7j to 8.3% (9.2%) for the optimistic (conservative) statistical estimate. To 
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conclude, a measurement of A with a precision of roughly 10% at the end of LHC 
running seems to be possible for our mAMSB benchmark scenario. 



5.4 Reconstruction of Couplings under the Assumption of a Wino LSP 

We have seen that the reconstruction of the coupling A is much more promising in 
SUSY scenarios with a wino like LSP than with a bino like LSP. Therefore, we will 
henceforth concentrate in this section on wino LSP scenarios. We will investigate 
the precision to which A can be measured as a function of the (left-handed) squark 
mass and the Xi mass. 

For that purpose, we perform a two-dimensional parameter scan. We have varied 
the pole masses of the left-handed first generation squarks (winos) between 400 GeV 
and 1210 GeV (100 GeV and 400 GeV). All other SUSY particles are decoupled. This 
part of parameter space has the chargino-squark-quark couplings equal to A = g to 
better than per mille precision. We employ a grid of ten times ten points equally 
spaced in the squark-wino mass plane. We assumed an integrated luminosity of 100 
fbT 1 at y/s = 14 TeV. For the relative error from the luminosity we use again 3%, 
whereas the errors from PDFs, NLO corrections, and the SUSY mass uncertainties 
were estimated in the following way. We calculated for a sub-grid with four times four 



parameter points the respective errors as it was done in Section [5]3] and extrapolated 
the errors for the other parameter points. 

The PDF uncertainties are described well by a linear function of the squark and 
Xi mass sum in which the error decreases with increasing mass. This behaviour can 
be understood by noting that the error of the gluon PDF decreases in the relevant 
region for increasing Bjorken-x [76]. For the NLO uncertainties we use a logarithmic 
fit as a function of the sum of the squark and Xi mass. Prospino employs the squark 
and Xi mass sum divided by two as its default factorization and renormalization 
scale. For the scale variation both scales are then varied between 1/2 and 2 times 
the default value, cf. Section [0| A logarithmic behaviour is expected, because the 
scale dependence enters via loops. 

The estimate of the error from SUSY mass uncertainties is more involved. In 
general, the precision to which the masses can be reconstructed depends on the 
nature of the SUSY particle, the mass spectrum and the model assumptions; see for 
example Refs. [4,84]. We here employ a simple approach. We assume that a squark 
mass of rriq L = 720 GeV can be reconstructed with an absolute error of Sm,q L = 10 
GeV. For a Xi with a mass of 100 GeV we assume the same error. This corresponds 



to our assumption made in Section |5.3|. We then estimated the error for other squark 



masses m~ £ via 



5mi L = 5m, L x < x M P * ± ||> , (5.2) 
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Figure 7: Lower cut for the jet-p^ (left-hand side) and upper cut for the jet invariant 
mass (right-hand side) as a function of the squark and Xi mass. The red region in the right 
figure corresponds to no cut on m(jetl). 



i.e. the relative error scales with the square root of the inverse squark pair production 
cross section, a(PP — > (? L (f L ), as expected if the errors are dominated by statistical 
uncertainties and if the masses are reconstructed from cascade decays of squarks. 

We also employ the same absolute error for the Xi mass, because the squark and 
Xi mass errors are always strongly correlated to each other if the squark decays into 
the Xi P ms a j e t [84]. Moreover, analyses such as the benchmark point SPSla [85] 
show that the lightest neutralino mass can be slightly better reconstructed than the 
squark masses [4, 84] . We therefore follow a conservative approach in this analysis. 
From the mass errors, we obtain errors for the total signal cross section following the 
procedure described in Section |5l| We found that a quadratic fit as a function of 
the squark mass describes very well the calculated errors as might be expected from 
eq. O). 

In order to find the right cuts for each SUSY scenario, we have employed a similar 
strategy to the one in Section |5~3| . We applied the first three cuts of Table ||. Starting 
from the trigger cut of 100 GeV, we searched numerically in steps of 10 GeV for a 
(lower) pt cut on the hardest jet that maximizes the significance S/ y/B, where 5* 
(B) corresponds to the number of signal (SM background events). Finally, we varied 
the (upper) cut on the invariant mass of the hardest jet in steps of 10 GeV in order 
to maximize the signal to SUSY background ratio. However, we only employ a cut 
on the invariant mass as long as the (conservative) significance is larger than five. 

The resulting cuts employed are shown in Fig. [7] for the jet-pr (which is equal 
to the cut on ^/ T ) and for the jet invariant mass. They are given as a function of 
the left-handed up-squark mass, rriu L , and the (wino-like) Xi LSP mass, m^o. Note 
that the left-handed down squark and left-handed up squark are quasi-degenerate, 
as are the xf an d Xv We can see that the lower cut on the jet-p^ increases with 
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Figure 8: Ratio of the number of signal events to the number of SUSY background 
events as a function of the squark and x!i mass. 



increasing squark mass and decreasing Xi mass as expected from the position of the 
Jacobian peak given in eq. ( |1 . 1|) . Recall that we usually obtain the best significance 



if we apply the jet-py cut close to the Jacobian peak; cf. Sections [4.2| and |572 . 

However, an exception is the quite large jet-p-r cut m the upper left corner of 
Fig. |7|. Due to the small squark-^ mass difference, the Jacobian peak lies below 
< 100 GeV, i. e below the trigger cut. For larger p^ values (compared to the peak) , the 
jet-pr (and j/ T ) distribution of the signal falls off less steep than the SM backgrounds. 
Therefore, for a small squark-%5 mass difference the best cut will lie above the 
Jacobian peak. 

Fig. |7| (right side) shows the upper cut on the invariant mass of the monojet. 
The red region corresponds to parameter points, where the conservative significance, 
S/ V7B, is smaller than five. In this case, we do not apply a cut on the invariant jet 
mass in order not to suppress the signal cross section. We observe that harder cuts 
on the invariant mass are usually applied for smaller squark and Xi masses, because 
the number of signal events is large and a hard cut will not suppress the signal below 
the discovery reach. Recall from Fig. f| and Fig. |6| that harder cuts on the invariant 
mass lead often to a better signal to SUSY background ratio. 

The importance of the invariant mass cut can be seen in Fig. ||, where the ratio 
of signal to SUSY background events (after cuts) as a function of the squark and Xi 
mass is shown. We observe that the ratio varies between 1.2 and 12. This good signal 
to SUSY background ratio justifies our approximation to neglect the (unknown) error 
from SUSY backgrounds as long as the SUSY background can be determined to a 
precision of a few 10%; see also the discussion in Section 



and Section [573. 

We also see some interesting structure in Fig. ||: For fixed squark mass, the sig- 
nal to SUSY background ratio first increases when we increase the xi mass and then 
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Figure 9: Fractional precision to which the XiQlQ coupling A can be reconstructed 
as function of the squark and x5 mass. The left (right) figure employs our optimistic 
(conservative) estimate for the SM background uncertainties. The solid and dashed black 
lines correspond to S/y/B (S/V7B) of 5a and 10<7, respectively. 



decreases again. For small xi mass, the SUSY background is dominated by Drell- 
Yan wino pair production. However, when the Xi mass is increased the respective 
cross section will decrease and squark pair production will take over to be the dom- 
inant background. At the same time, as long as the Xi mass increases moderately 
(relative to the squark mass) the signal cross section decreases less fast than the 
SUSY background leading to a better signal to SUSY background ratio. However, 
when we come to the region where m^o w we obtain a reduced signal to SUSY 
background ratio. The signal to SUSY background ratio can be further increased 
with a cut on the jet invariant mass. As can be seen in Fig. || this can result in a 
signal to SUSY background ratio larger than ten! 

The fractional precision that can be achieved can be seen in Fig. [| where we 
show the relative error of the coupling A as a function of the left-handed up-squark 
mass and the Xi LSP mass. Fractional precisions between 0.12 and 0.24 are possible 
throughout parameter space for optimistic assumptions, or 0.12-0.44 for conservative 
ones. The most precise measurements are for lighter sparticles, where signal event 
numbers are higher. 

The significance of the SUSY monojet signal is also given in Fig. |J] by the solid 
(dashed) black line corresponding to the 5a (10a) region. We see a strong dependence 
on whether we assume optimistic or conservative error estimates. Whether or not 
the signal can be seen over statistical fluctuations of the background to the 5a level, 
the coupling A may be bounded. 

To conclude, the LHC will be able to test the XiQlQ coupling relation to a 
precision of O(10%) in large regions of the wino LSP parameter space as long as the 
squark masses are < 1 TeV. 
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6. Summary and Conclusion 

If signatures of new physics beyond the SM are discovered that are consistent with 
SUSY, a program of measurements of the interactions of new states would be useful 
to hypothesis test the SUSY interpretation and cross-correlate with other measure- 
ments. As a first step, we have investigated the feasibility of SUSY monojet produc- 
tion at the LHC via qg — > qxi — > qXiXi f° r measuring the Xa.QQ coupling A, under 
the assumption that the Xi is the LSP. 

For a bino LSP, supersymmetry predicts that A is proportional to g', the U(l)y 
gauge coupling. For a wino LSP however, supersymmetry predicts A = g. When 
added to other empirical information on the field content of the neutralino, the 
value of the coupling constitutes a test of supersymmetry. In general, the MSSM 
Xi contains admixtures of the wino, the bino and the higgsino, resulting in A being 
a mixture of g, g' and 0, respectively. Without additional information on the field 
content coming from other measurements, a measurement of the coupling of the LSP 
to a quark and a squark (in the framework of the MSSM) yields valuable information 
regarding the neutralino field content. As such, it could be used in the calculation of 
its thermal relic density in the universe, or for dark matter direct detection rates, for 
example. Our signal would only be measured significantly above backgrounds after 
data have been amassed for many years at the LHC. Therefore it is reasonable to 
expect that information on the mass spectrum would already be known, e.g. squark 
and neutralino mass scales, thereby constraining the SUSY backgrounds. 

At this stage, it would be useful to extract, as far as possible, SUSY backgrounds 
from other data samples. However, without data to guide us, the number of possi- 
bilities is very large. The most difficult part of estimating SUSY backgrounds would 
likely be the extraction of the relevant branching ratio for the q — > qx\ decay. Nev- 
ertheless, placing an upper bound of 1 on the branching ratio for an observed SUSY 
monojet signal would yield a lower bound on A. If the LSP is a generic admixture 
of bino-wino-higgsino eigenstates, one must also take into account the fact that the 
lightest chargino may couple differently to the lightest neutralino. However, one 
would then generically expect the Xi to not be quasi-mass degenerate with the x^, 
and so its decays ought to be visible and therefore reducible by our cuts. On the 
other hand, it may be that a simple model of supersymmetry breaking with only a 
few parameters is selected by fitting all available data. In this case, the branching 
ratio and therefore the total rate of our monojet signal would be predicted and would 
constitute an additional hypothesis test on that model. 

Measuring A is not easy at the LHC, but we have shown that a fractional accuracy 
approaching 10% is feasible with a large data set, providing a non-trivial and useful 
empirical constraint on the MSSM. As such, it will be an important element of the 
proposed program of precision measurements of new particle interactions to support 
the SUSY hypothesis. 
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A. The Jacobian Peak in the Transverse Momentum Distri- 
bution 

The Jacobian peak is a well-known feature of the transverse momentum distribution 
of the electron in the process A + B — y W ± + X — y e ± + v + X, where A and B 
are the initial state hadrons. The resulting peak at pr — \m w is a consequence 
of the Jacobian that arises from changing kinematic variables from cos# (where 6 
is the center-of-mass scattering angle) to pr- 10 In this paper, we have focused on 
monojets that arise from qxi production, where q — > qx\, and the quark is observed 
as a hadronic jet. The pt distribution of the quark jet also exhibits a Jacobian peak. 
In this appendix, we derive an approximate expression for the location of the peak 
in the transverse momentum distribution of the jet. 

Consider the 2 — y 3 scattering process, which schematically is of the form: 

a + b ->■ c + 3, followed by c -» 1 + 2 , (A.l) 

where the decaying particle c is spinless. Since the particles a, b and 1 represent light 
quarks or gluons, we shall set their masses to zero, m a = = mi = 0. We denote 
the mass of particle c (identified as the q) to be m c = M, and the masses of particles 
2 and 3 (which are identified with Xi) t° be rri2 = = m. 

If the particle c is on-shell, then the corresponding matrix element for the 2 — y 2 
process, a + b — )>c + 3isof the form 

M(a + b ->■ c + 3) = d(s,t) , (A.2) 
10 For a pedagogical treatment, see Ref. [86]. 
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where Ci(s,t) is a dimensionless function of s = (p a + Pb) 2 , t = {p a — p c ) 2 and the 
particle masses. The kinematical limits of t are: 

_I [ s -M 2 -m 2 + \ 1/2 (s, M 2 , m 2 )] <t<-\[s-M 2 -m 2 - X 1 ' 2 ^, M 2 , m 2 )] , 

(A.3) 

where 

X(a, b, c) = a 2 + b 2 + c 2 - 2a6 - 2ac - 26c (A.4) 

is the well-known triangle function of relativistic kinematics. The squared matrix 
element for the decay of particle c (which is either q^ or qn), summed over final 
spins, is given by 

\M(c ^ 1 + 2)\ 2 = C 2 (M 2 - m 2 ) , (A.5) 

where C 2 is a dimensionless (real positive) constant that will eventually cancel out in 
our computation. Using, eq. ( |A.5| ) it follows that the total width of particle c times 
the branching ratio is given by 

C 2 M ( m 2 \ 2 / \ \ 

Br = -^-(i--r 2 ) , (A.6) 



16tt V M 2 f 

where the branching ratio B = B(c — > 1 + 2). 

To set up our computation, we work in the center-of-mass system. Then, the 
four-vectors of the initial states and the observed final state (particle 1) are: 

p tt = |vS(l; 0,0,1), (A.7) 
P6 = !VS(1; 0,0,-1), (A.8) 
pi = £i(l; sin^, 0, cos0) , (A.9) 

where 9 is the scattering angle in the center-of-mass frame. Following Ref. [87], we 
define four Lorentz-invariant quantities, 

t 1 = (p a -p 1 f = -^- s E x {l-cos9), (A. 10) 

t 2 =t = (pi + P2 -p a ) 2 , (A. 11) 

Sl =p 2 c = (p 1+ p 2 ) 2 , (A.12) 
s 2 ee (p a +p b - pi) 2 = s-2y/sEx. (A.13) 

We denote the three-body phase space integral by 

R 3 (s)= I \[ < ^5^\p a +p b -p 1 -p 2 -p z )5{^- s -E l -E 2 -E,). (A.14) 
J i=i % 

The key formula that we need is given by eq. V-7.8 of Ref. [87], 

Q{-G(s, h, 82, m 2 a , m 2 b , m{)} 



^ 3 71 rvf nl 4. 2 2 2\ 



ds 2 dtidsi 8A 1 / 2 (s, m\, m^)A 1 / 2 (s, s 2 , m\) 



xQ{—G(si,s 2 ,s,m 2 ,m 1 ,m 3 )} / dcj) , (A. 15) 
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where G is the basic four-particle kinematic function first introduced in Ref. [88], 



I / 2u x + u — v u + w — y\ 

G(x, y, z, u, v, w) = — - det I x + u — v 2x x — z + w J , (A.16) 

\u + w — y x — z + w 2w / 

and A is the triangle function defined in eq. (|A.4j) . Expanding out the determinant 
yields the unwieldy expression, 11 

G(x, y, z, u, v, w) = xy(x + y) + zu(z + u) + vw(v + w) + x(zw + uv) + y{zv + uw) 
—xy{z + u + v + w) — zu{x + y + v + w) — vw(x + y + z + u) . (A. 17) 

Finally, is the so-called helicity angle [87], which is most conveniently defined in 
a reference frame where p 2 + p 3 = p a + p b — p x = 0. In this reference frame, is 
identified as the azimuthal angle between the production plane spanned by p b and 
p 1 and the plane spanned by p x and p 3 , with p x as the axis. The angle 4> can be 
re-expressed in terms of the Lorentz-invariant variables S\, s^, t\ and tz, as exhibited 
in eq. V-8.8 of Ref. [87]. In particular, it will be convenient to express £2 in terms of 
si, s 2 , ti and cos0 following eq. V-8.9 of Ref. [87], 



ml + m\- t 2 



D — 2 [G(s,t\, S2, ml, ml, m\)G{si, S2,s, m\,m\, m^)] 1 ^ 2 cos0 

\(s,s 2 ,ml) 

(A.18) 



where 




s + s 2 — m\ s — Si + ml 

2s 2 s 2 - m\ + m\ | . (A.19) 

s 2 - ti + ml 

Note that the phase space distribution in the helicity angle is uniform, as the inte- 
gration over <j) in eq. ( |A.15| ) is trivial. However, because the matrix element given 
in eq. ( |A.2| ) depends on t = t 2 , the calculation of the partonic cross section for 
a + b — > 1 + 2 + 3 will require a nontrivial integration over 0. 

The step functions in eq. ( |A.15 ) determine the kinematical ranges of the pa- 
rameters Si, S2 and t\. Taking m a = m b = mi = and m 2 = m 3 = m, it follows 
that: 



ds2dtidsi 8s(s 



r2n 

r e{-G(s, h, s 2 , 0, 0, 0)}9{-G(ai, s 2 , s, m 2 , 0, m 2 )} , 

-82) JO 



(A.20) 



n Eq. (|A.17|) , which was first defined in Ref. [88], is also given in eq. IV-5.23 of Ref. [87]. We 
have noted a typographical error in the latter; in the second line of eq. IV-5.23, the first term yzw 
should read yzv. 
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The differential cross-section is given by: 

f 



da 



647T 5 S 



dR 3 (s)\M(a + b -+ l + 2 + 3)|' 



(A.21) 



where the squared matrix element is suitably averaged over initial spins and summed 
over final spins. The dominant contribution to a + 6 — > 1 + 2 + 3 takes place via 
a + b — > c + 3, where c is produced approximately on-shell and subsequently decays 
via c — > 1 + 2. In particular, since c is a spin-zero particle, 



\M{a + b^ l + 2 + 3)| 



2 _ \M(a + b^ c + 3)\ 2 \M(c^ 1 + 2)| 2 
(sj - M 2 ) 2 + M 2 r 2 



(A.22) 



We now use eqs. ( |A.2| ) and ( |A.5| ) and employ the narrow width approximation 

1 7T 



(si -M 2 ) 2 + M 2 r 2 



Mr 



<5(si - M ) . 



(A.23) 



Hence, it follows that: 
da 



B 



ds 2 dt\dsi 32ti 2 ^s 2 {s — s 2 



-8{ Sl -M 2 )e{-G(s, t 1>S2 , 0,0,0)} 



xO{-G(s 1 ,s 2 ,s,m,0,m 2 )} / \C 1 (s,t 2 



2n 



(A.24) 



where ^2 should be expressed in terms of si, S2, ti and cos0 using eq. ( |A.18| ) before 
performing the integration over 0, and 



m 
M 2 



(A.25) 



Assuming that G(si, S2, s, m 2 , 0, m 2 ) < 0, we can immediately use the ^-function 
to integrate over Si. Using eq. ( |A.17|) , we obtain: 



G(si, s 2 , s, m 2 , 0, m 2 ) = s 2 s 2 - s 1 s 2 (s - s 2 + 2m 2 ) + m 2 s(s - s 2 ) + s 2 m 4 

= s 2 (s 1 -st)(s 1 -s^), (A.26) 



where S2 is strictly non-negative and 



o 1 —m 2 + |(s - s 2 ) 
That is, we require that: 



1 + Wl 



4m 2 
s 2 



(A.27) 



sr < M 2 < sf , (A.28) 
otherwise, s\ = M 2 can never be satisfied when G(s\, s 2 , s,m 2 , 0, m 2 ) < 0. Note that 



eq. ( |A.28j ) yields upper and lower limits for s 2 . One can then use eq. (|A.13 ) to obtain 
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upper and lower limits for E\. These limits correspond to the roots of the quadratic 
equation, 

4^M 2 Ef - 2(M 2 - m 2 )(s + M 2 — m 2 )E x + a/s(M 2 - m 2 ) 2 = . (A.29) 
These roots can be expressed as: 12 



Ef = -4= \s + M 2 - m 2 ± A 1/2 (s, M 
4v/i 



2 ,m 2 ), . 



(A.30) 



where £ is defined in eq. ( |A.25| ). Likewise, employing eq. ( |A.13| ), we define 

sf = s- 2^Ef . (A.31) 
The range of t\ is determined from the inequality: 

G{s, t u s 2 , 0, 0, 0) ee sti{s + h - s 2 ) < , (A.32) 



where we have used eq. ([A. 171) to evaluate the G-function. That is, as s 2 ranges over 
s 2 < s 2 < s%, 

s 2 -s<tx<0. (A.33) 



Assuming that s 2 < s 2 < , the integration of eq. ( |A.24| ) over sx is immediate, 
and we obtain: 

t^t = QO 2C f re{-G( s ,t 1)S2) o,o,o)} / id^A + ^cos^)! 2 ^, 

(A.34) 

where eq. (|A.18|) has been used to write t 2 = A 1 + A 2 cos <p. Using eqs. ( |A.26|) and 
( |A.32 ) with s\ — M 2 , m a = nib = rni = and m 2 = = m, the coefficients A\ 
and A 2 are given by 



A x = m - 



, s 2 (s - s 2 + ti)(M 2 - m 2 ) - stx(s - s 2 - M 2 + m 2 ) 



(A.35) 



{s-s 2 y 



.4, 



2 [sti(s - s 2 + ti)(M 4 s 2 - M 2 s 2 (s -s 2 + 2m 2 ) + m 2 s(s - s 2 ) + s 2 m A )\ 



[s - s 2 



(A.36) 



We now introduce the transverse momentum, p? of particle 1, which is defined 
by Pt = Ex sin 6. Note that 



stx(s + tx — s 5 



2 2 
-3 p T 



(A.37) 



Note that eq. (A.3C) is equivalent to Ef = ^(E c ±p c ), where E c andp c are the center-of-mass 



energy and momentum of the decaying particle c. 
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which is strictly non-positive as required by eq. ([A. 32 ). In particular, 

cos0 = ±4/l- §r (A.38) 



where the ± indicates that 9 and it — 6 correspond to the same value of px- Thus, 
eq. flATTOl) yields 

(A.39) 



One can now perform a change of variables from {t\ , s 2 } to \p\ , E{\. Comput- 
ing the Jacobian of the transformation, it follows that: 

, , sdp%dEi ,. lnS 

dt 1 ds 2 = P J 1 . (A.40) 

V E{ - Pt 

The limits of the kinematic variables px and Ei are given by: 

< p T < Ei , £f < Ei < E+ , (A.41) 

where the range of pt follows from | cos6\ < 1 and Ef is defined in eq. ( |A.30| ). Since 
we aim to compute da/dpr, it is more useful to interchange the order of integration. 
Thus, equivalent to eq. ( |A.41| ) is: 

for 0<p T <Ei , Ei < E 1 < Ef , (A.42) 

for Ei <p T < Ef , Pt<E x < ■ (A.43) 

Combining eqs. (|A.34| ) and (|A.4Q ), and adding the contributions from the two pos- 
sible values of ti that yield the same value of pt [cf. eq. ( |A.39| )1, one obtains: 

U ^ d<P V|d( S ,4 j) +A 2 cos0)| 2 , (A.44) 



dp\dEi QATi 2 is^ 2 Ei^/El-p 2 T 



where is defined by eq. (|A.35|) , and the superscript indicates which sign is used in 
eq. ( |A.39|) to express ti in terms of Ei and p T . In contrast, A 2 [defined in eq. ( |A.36|) 1 
does not depend on the sign choice in eq. (|A.39|) as a consequence of eq. (|A.37|) . 

We now integrate over Ei, employing the limits of integration given in eqs. ( |A.42|) 
and ( |A.43j ). Writing dp\ = 2p?dpT, we arrive at our final result, 13 

dp T 32vr 2 es 3 / 2 J Emin Ei^/E\-p\ J 



where the upper and lower limits of integration are given by E miiyi = E± and 

E m ,„=r' for0£ ^ £ ^' + (A.46) 
for Ei <p T < Ef . 



There is no singularity in the limit of £ — > since in this limit, E ± -t 0. 



- 35 - 



As a warmup, we shall ignore the details of the scattering matrix element for the 
process a + b — )• c + 3 by putting G\ = 1. In this case, the integrals in eq. (|A.45|) are 
elementary, and the end result is: 



da 



B 



dp T 8tt£s 3 / 2 



tan 



Pt 



Q(Ei -p T ) tan" 1 



"12 



Pt 



Pt 



where < pt < E±, and the step function 6 is defined as usual, 
Q{E^-p T ) = 



for < p T < E^ , 
for El <p T < E+ . 

It is convenient to introduce dimensionless variables, 



(A.47) 



(A.48) 



w = 



2E 1 



x = 



2pi 



y 



M 2 
s 



m 

w 



(A.49) 



The kinematics of the scattering process requires that yfs > M + m, which is equiv- 
alent to the condition, 

y/y{l + Vz)<\. (A.50) 
The range w is given by w~ < w < w + , where 
,± 



w~ 



|(1 -z) [l + y{l-z)±X 1 / 2 (l,y,yz 



Hence, the range of x is 



< x < w + < 1 



(A.51) 



(A.52) 



As an example, take y = 0.5 and z = 0.1, which is consistent with the inequality 
given in eq. ( |A.50|) . Eq. ( |A.52[ ) then implies that < x < 0.79657. The transverse 
momentum distribution, plotted in Fig. |10| exhibits a striking Jacobian peak located 
at x = 0.50843, which corresponds to 



T ) peak 



El- 



(A.53) 



The origin of the Jacobian peak is a consequence of the change in kinematic 
variables given in eq. QA.40 ), and is rather insensitive to the form of the matrix 
element. To illustrate this point, we have numerically evaluated eq. ( A.45Q , where 
the tree-level form for C\ for gq — > qnXi is employed [11,31,52] 



|Ci(M)| 2 = iV 



s + t- M 2 M 2 (m 2 -t) sm 2 + (m 2 - t)(M 2 - m 2 ) 



(A.54) 



2s (M 2 -t) 2 s(M 2 -t) 

where is an overall dimensionless normalization factor that depends on the relevant 
couplings. In terms of the dimensionless variables introduced in eq. QA.49 ), 



da Bx 

dx 647r 2 s(l — z) 



/ r=2 2 / d^22 C 1 {s,A ( f ) +A 2 

J^in W VW 2 ~ X 2 J ^ 



cos0)| 2 , (A.55) 
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Figure 10: Unnormalized pt distributions for a + b — > c + 3 , c — > 1 + 2, assuming that 
the matrix element for a + b — )• c + 3 is constant (dashed curve) or is given by eq. ( |A .54| ) 
(solid curve). The rescaled transverse momentum is defined by x = 2pT/yfs and can take 
on values in the range < x < x max , where x max = 3(1 — z) [l + y(l — z) + A 1 / 2 (l, y, yz)] . 
The masses of particles c and d are fixed by y = M 2 / s = 0.5 and z = m 2 /M 2 = 0.1. 
To facilitate the comparison of the two pt distributions, the relative normalization of the 
two curves has been fixed such that the height of the peaks of the distributions coincide. 



The location of the peak at x = 0.50843, corresponding to eq. ( A. 53 ), is the same for both 
curves. 



where 



w 



UK 



W 



and the coefficients and A2 are given by: 



(±) 



\s\y{l + z)-l± 



y/w 2 



for < x < w , 
for w~ < x < w~ 



+ - z)(w - 2)] 



An 



sxy 



1/2 



XV, 



[w(l — w — z) — 1/(1 — w)(l — z) 2 ~\ 



21 1/2 



(A.56) 



(A.57) 
(A.58) 



The resulting unnormalized pt distribution is exhibited in Fig. [10]. Note that 
the shape of the pt distribution is dominated by the explicit kinematic factors that 
appear in eq. ( |A.45|) , and depends quite weakly on the actual form of the squared- 
matrix element given in eq. ( |A.54|) . Moreover, the location of the peak in the pt 
distribution is unchanged and given by eq. ( |A.53 ), as a consequence of structure of 
the kinematic limits given in eqs. ( |A.42j ) and ( |A.43| ). 

In the above analysis, the location of the Jacobian peak given in eq. ( |A.53| ) 
depends on the partonic center-of-mass energy The differential cross section for 



-37- 



the hadronic scattering process, A + B — > c + 3 + X — > 1 + 2 + 3 + X, is obtained by 
convoluting the pr distribution of the partonic subprocess, a + b — > c + 3— > 1 + 2 + 3, 
with the product of the parton distribution functions fai x iyQ 2 )f^( x 2,Q 2 ), where 
the total center-of-mass squared-energy S is related to the partonic center-of-mass 
energy via s = X1X2S, and Q is the factorization scale. In the convolution, partonic 
center-of-mass energies close to the energy threshold for the partonic process provide 
the dominant contribution to the production of the final state. In this case, one 
can derive an approximate formula for the location of the Jacobian peak that does 
not depend on the partonic center of mass energy. The threshold for a + b — > c + 3 
corresponds to the point at which 

A(s, M 2 , m 2 ) — (s + M 2 — m 2 ) 2 - AsM 2 = . (A.59) 

At this point s + M 2 — m 2 = 2My/s (or equivalently, y/s = M + m), in which case 

, M 2 -m 2 , K , 

E ^ = E f = ^M-- < A ' 60 > 



Of course, the cross-section given in eq. (|A.45| ) vanishes exactly at threshold where 
E± = . However, if we are close to threshold, then eq. ( |A.60|) still provides a 
decent approximation to E± , in which case the location of the Jacobian peak is: 

Mpeak = #r - 2M = h£M , (A.61) 

which is independent of the partonic center-of-mass energy. 

In this paper, we have numerically computed the transverse momentum distribu- 
tion of the hadronic scattering process, taking into account the partonic scattering 
process at all allowed values of the partonic center-of-mass energy. In particular, 
as the partonic center-of-mass energy is increased above the threshold energy for 
a + b — > c + 3, the location of the peak of the partonic transverse momentum dis- 
tribution, E^ [cf. eq. ( |A.30| )1 decreases relative to the estimate given in eq. (|A.61| ). 
Thus, we expect the actual peak in the transverse momentum distribution of the 
hadronic scattering process (or equivalently in the missing transverse energy distri- 
bution) to be somewhat less than the result of eq. ( |A.61 ). This is indeed the case in 
the j/ T distributions that we exhibit in this paper. 

Note that in the approximation that the transverse momentum of particle c is due 
entirely from the hard scattering process (i.e. the transverse momentum of the initial 
partons and the spectators are neglected), the distribution of the missing transverse 
energy (i.e. particles 2 and 3 of the hard scattering process) should precisely match 
that of the transverse momentum of the monojet (i.e. particle 1 of the hard scat- 
tering process). Of course, the effects of spectators, initial and final state radiation, 
fragmentation of final state partons, jet mismeasurements and detector effects will 
tend to reduce the sharpness of the peak in the jj T distributions as compared to that 
of Fig. 0. 
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